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Abstract.  The  characteristics  of  the  hypersonic  flow  around  an  elliptical- 
cross  section  cone  and  the  computational  code  to  perform  large-eddy  sim- 
ulations of  this  flow  are  described.  Preliminary  results  of  the  Lagrangian 
implementation  used  to  compute  the  subgrid-scale  terms  are  presented. 


1.  Introduction 

Many  aspects  of  transitional  and  turbulent  flows  are  not  fully  understood. 
This  is  especially  true  in  the  hypersonic  regime,  where  examples  of  un- 
resolved issues  include  the  effect  of  freestream  disturbances  and  three- 
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dimensionality.  In  the  absence  of  detailed  experimental  or  computational 
databases  to  better  understand  these  physical  phenomena,  we  are  left  with 
excessive  design  conservatism  and  unrefined  conceptual  designs. 

When  investigating  these  phenomena  via  CFD,  direct  numerical  simu- 
lations (DNS)  are  not  affordable.  Turbulence  models  provide  a wide  range 
of  accuracy  in  predicting  turbulent  flows  of  engineering  interest.  Depend- 
ing on  the  level  of  detail  required,  one  may  choose  the  Reynolds-averaged 
Navier-Stokes  (RANS)  models  or  the  state  of  the  art  subgrid  scale  (SGS) 
models  in  a large-eddy  simulation  (LES)  to  obtain  a more  refined  predic- 
tion. The  study  of  fundamental  physical  phenomena  must  be  done  using 
the  finest  possible  prediction,  namely  LES.  However,  one  must  keep  in  mind 
that  a key  feature  of  the  prediction  is  the  validation  with  experiments. 

Using  the  most  recent  laser  and  camera  technologies,  Huntley1  and 
Huntley  et  al.2  present  the  first  detailed  flow  visualization  of  natural  tran- 
sition on  an  elliptical  cross-section  cone  at  Mach  8.  Mean  flow  features  and 
details  about  the  unstable  modes  in  the  boundary  layer  for  the  same  con- 
figuration are  given  by  Kimmel  et  al .3’4  and  Poggie  et  al.,5  respectively.  Be- 
cause this  flow  is  being  extensively  documented  experimentally  and  because 
the  geometric  configuration  resembles  that  of  the  forebody  of  a hypersonic 
vehicle,  this  database  is  ideal  to  test  the  state  of  the  art  SGS  models  and 
LES  methodology  for  hypersonic  flows. 

The  present  work  is  an  ongoing  effort  to  provide  detailed  flow  simula- 
tions of  unsteady,  hypersonic,  transitional  or  turbulent  flows.  In  Martin  et 
al.  (2000a)  we  develop  and  test  SGS  models  for  compressible  LES  using  the 
apriori  test  in  compressible,  isotropic  turbulent  flow.  In  Martin  et  al.  (2001) 
we  validate  the  LES  methodology  against  the  DNS  of  supersonic  boundary 
layers,  and  in  Martin  et  al.  (2000b)  we  extend  the  LES  code  to  general- 
ized curvilinear  coordinates  and  validate  the  implementation  in  supersonic 
turbulent  boundary  layer  flow.  In  this  paper,  we  present  preliminary  LES 
results  of  the  flow  around  a section  of  an  elliptical  cross-section  cone  away 
from  the  tip  of  the  cone.  The  flow  conditions,  simulation  procedure,  prelim- 
inary flow  assessments,  and  future  work  are  given  in  the  following  sections. 

2.  Geometry  and  flow  conditions 

Following  the  experimental  configuration  of  Huntley  et  al.2  the  lifting  body 

geometry  is  a sharp-nosed,  elliptical  cross-section  cone.  The  nominal  radius 

of  the  experimental  model  is  less  than  200  pm.  In  our  simulations,  we  use  an 
ellipsoidal  nose  with  an  80  pm  diameter  of  nose  ellipsoid  measured  on  the 
major  axis.  The  afterbody  is  an  elliptical  cross-section  cone  of  4:1  aspect 
ratio,  17.5  degrees  half-angle  in  the  major  axis,  and  0.242  m in  length, 
resulting  in  base  dimensions  of  0.152  m across  the  major  axis  and  0.038  m 
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across  the  minor  axis. 

The  freestream  flow  conditions  are  Re oo  = 14  x 106  /m,  M0 0 = 8, 
Poo  = 0.5,  and  = 58  K.  For  these  conditions,  the  boundary  layer  is  fully 
turbulent  at  position  x=  1 7.5  cm  from  the  nose.1  The  wall-temperature 
condition  is  prescribed  to  450  K,  which  is  nearly  adiabatic.  As  in  the  ex- 
periments, air  is  the  working  fluid. 

The  Mach  8 flow  around  the  4:1  elliptical  cross-section  cone  at  zero  angle 
of  attack  is  highly  compressed  behind  the  shock.  The  difference  in  shock 
strength  between  the  major  and  minor  axis  causes  a higher  compression 
at  the  leading  edge  producing  a cross  flow  from  the  leading  edge  to  the 
centerline  of  the  cone.  This  is  illustrated  in  Fig.  la.  At  the  centerline  the 
cross-flow  velocities  are  zero  and  mass  conservation  induces  a bulge,  see 
Fig.  lb,  where  the  boundary  layer  is  twice  as  thick  as  the  boundary  layer 
in  the  off-center  region.1  Experiments  show  that  transition  occurs  first  on 
the  centerline.1 

3.  Grid  resolution 

We  estimate  the  resolution  requirements  by  considering  the  dimensions  of 
the  turbulence  structure  and  the  data  provided  by  the  experiments.  From 
the  experimental  data1,  the  estimated  wall  unit  is  7.5  x 10-5  m and  the 
boundary  layer  along  the  centerline  is  laminar  at  x = 9 cm,  late-transitional 
at  11.4  cm,  and  fully  turbulent  at  17.5  cm  from  the  nose.  See  Figure  2.  The 
turbulent  boundary  layer  at  the  center  line  is  5 to  6 mm  thick.  Note  that 
the  shock-standoff  distance  is  6S  = 1.2 9C,  where  9S  and  9C  are  the  shock 
and  cone  angles,  respectively.  The  turbulence  structures  on  a flat  plate 
boundary  layer  extend  about  100  wall  units  in  the  spanwise  direction  and 
a few  boundary  layer  thicknesses  in  the  streamwise  direction. 

In  the  spanwise  direction  we  use  a uniform  grid  spacing  of  22  wall  units. 
This  estimate  is  based  on  the  grid  resolution  used  in  previous  LES  of  a 
supersonic  boundary  layer. 7 For  the  conditions  chosen,  the  boundary  layer 
at  the  leadingedge  is  laminar.1  Since  the  flow  is  supersonic  and  the  cone  is 
at  zero  angle  of  attack,  the  flow  around  the  top  of  the  cone  is  not  affected 
by  the  flow  on  the  bottom.  Thus  only  the  top  of  the  cone  is  simulated.  In 
the  streamwise  direction,  the  flow  is  laminar  at  x = 9 cm  from  the  nose. 
Thus,  the  resolution  requirements  up  to  9 cm  from  the  nose  are  given  by 
the  grid  convergence  studies  of  the  laminar  flow  at  the  nose.  From  x = 9 
cm  to  the  base  of  the  cone  we  require  that  the  maximum  grid  spacing 
on  the  surface  is  33  wall  units.  In  the  wall-normal  direction  the  grid  is 
exponentially  stretched,  we  require  0.3  and  S'"  wall  units  for  the  minimum 
and  maximum  grid  spacings  within  the  boundary  layer,7  where  <j+  = S/zT 
is  about  70  for  the  turbulence  case.  The  computational  domain  is  large 
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enough  to  include  the  standoff  shock  wave. 

4.  Simulation  procedure 

For  the  LES,  we  use  a third-order  accurate  WENO9  to  compute  the  convec- 
tive fluxes.  This  scheme  has  low  dissipation  properties,  and  was  designed 
to  perform  DNS  and  LES  of  compressible  flows.  The  time  advancement 
technique  is  based  on  the  DPLU  relaxation  method  of  Candler  et.  a/10  and 
was  extended  to  second-order  accuracy  by  Olejniczak  and  Candler.11  The 
viscous  fluxes  are  evaluated  using  fourth-order  central  differences.  Finally, 
the  transformation  metrics  are  evaluated  using  fourth-order  central  differ- 
ences so  that  the  inaccuracy  of  the  numerical  evaluation  of  the  metrics 
coefficients  is  less  than  the  inaccuracy  of  the  convected  fluxes. 

The  initial  condition  will  be  a superposition  of  laminar  flow  and  a pre- 
scribed freestream  energy  disturbance  spectrum.  To  generate  the  laminar 
solution  we  use  a finite  volume  code,12  where  we  only  compute  90  degrees 
of  the  cone  geometry  and  use  bilateral  symmetry  to  reflect  the  resulting 
solution  across  the  centerline  and  generate  the  full  180  degrees.  We  then 
interpolate  the  laminar  solution  from  the  finite-volume  cell  centers  to  the 
finite-difference  grid  points  using  tri-linear  interpolation. 

The  required  SGS  terms  and  the  model  representations  are  given  in 
Martin  et  al.  (2000a,  2000b).  To  evaluate  the  model  coefficients  we  use  the 
Lagrangian-averaging  operation,  where  the  averaging  is  performed  along  a 
fluid  particle  pathline.  A full  description  of  the  Lagrangian-average  proce- 
dure can  be  found  in  Meneveau  et  al.  (1996). 

5.  Preliminary  flow  assessments  and  future  work 

In  this  section  we  present  a brief  progress  report  in  performing  the  LES. 
A portion  of  the  cone  has  been  initialized  using  the  interpolated  laminar 
solution  to  test  the  Lagrangian  implementation  of  the  SGS  models.  To 
minimize  the  complexity  of  this  test,  the  disturbances  introduced  by  the 
tri-linear  interpolation  are  used  as  initial  disturbances.  Figures  3 through  5 
show  contours  of  the  SGS  terms  at  the  exit  plane  of  the  cone  on  spanwise 
wall-normal  planes.  These  figures  include  a quarter  of  the  computational 
domain  (centered  about  the  centerline)  in  the  spanwise  direction,  and  about 
ten  boundary  layer  thickness  in  the  wall  normal  direction,  the  boundary 
layer  thickness  along  the  centerline  is  about  2 mm,  three  times  smaller 
than  the  turbulent  boundary  layer.  Figures  3 through  5 show  that  the  SGS 
terms  are  dominant  in  the  bulge  region.  This  would  not  be  the  case  if  the 
model  coefficients  were  calculated  using  the  ensemble  average  procedure. 

The  magnitude  of  the  initial  disturbances  is  very  small  and  do  not  grow 
significantly  in  time.  Future  work  includes  imposing  a prescribed  distur- 
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bance  energy  spectrum  and  running  the  LES  to  gather  sufficient  statistical 
data  to  assess  the  effect  of  freestream  disturbance  and  study  the  transition 
and  turbulence  phenomena. 

Computer  time  was  provided  by  NASA  Ames  Research  Center.  This 
work  was  supported  in  part  by  the  Air  Force  Office  of  Scientific  Research 
under  grant  AF/F49620-01-1-0060. 
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Figure  1.  (a)  Pressure  and  (b)  density  contours  in  the  exit  plane  of  the  cone  geometry 

for  the  laminar  solution. 
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Laminar  at  9 cm 
Late  transitional  at  11.4  cm 


Fully  turbulent  at  17.5  cm 


24.15  cm 


3.81  cm 


Shock  stand-off  distance 
0.314  cm 


Shock  stand-off  distance 

1.62  cm  = 3.24<5  (using  5 = 5 mm  and  9S  = 1.20c) 


Figure  2.  Model  dimensions  and  key  flow  features  required  to  determine  the  grid  reso- 
lution and  computational  domain  size. 


Figure  3.  SGS  stress  contours  on  the  exit  plane,  the  flow  is  into  the  page. 
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Figure  4-  SGS  heat  flux  contours  on  the  exit  plane,  the  flow  is  into  the  page. 


Figure  5.  Contours  of  SGS  turbulent  kinetic  energy  on  the  exit  plane,  the  flow  is  into 
the  page. 


